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Abstract 

We study the angular power spectrum estimate in order to search for large scale 
anisotropies in the arrival directions distribution of the highest-energy cosmic rays. 
We show that this estimate can be performed even in the case of partial sky coverage 
and validated over the full sky under the assumption that the observed fluctuations 
are statistically spatial stationary. If this hypothesis - which can be tested directly 
on the data - is not satisfied, it would prove, of course, that the cosmic ray sky is non 
isotropic but also that the power spectrum is not an appropriate tool to represent 
its anisotropies, whatever the sky coverage available. We apply the method to sim- 
ulations of the Pierre Auger Observatory, reconstructing an input power spectrum 
with the Southern site only and with both Northern and Southern ones. Finally, 
we show the improvement that a full-sky observatory brings to test an isotropic 
distribution, and we discuss the sensitivity of the Pierre Auger Observatory to large 
scale anisotropies. 



1 Introduction 



The origin of the highest energy cosmic rays is a theoretical challenge of mod- 
ern astrophysics, and is subject of much experimental efforts. Above 10 20 eV, 
the current data are too scarce for one to make any definitive statement about 
the existence or the lack of the GZK cutoff, as well as a statistically meaning- 
ful information about the arrival direction distribution. Whereas the AGASA 
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experiment ( Takeda et al. 2003h is over since January 2004, a new generation 



of experiments especially dedicated to the highest energies is emerging, and 
the first of them is the Pierre Auger Obs ervatory currently under construc- 
tion ( The Pierre Auger collaboration 19 95). For a recent review on the state o f 



the art of the highest-energy cosmic rays, we refer the reader to ((Cronin 20041 ) 
for instance. 

The distribution of the arrival directions is certainly one of the most crucial 
observable in order to yield some evidences about the sources of the highest- 
energy cosmic rays ijlsola et al. 2002USigl et al. 200.4 Isid et al. 2004h . When 



trying to point out large scale anisotropics, one is naturally led to work with 
the angular power spectrum of the arrival direction distribution. The detec- 
tion of large scale anisotropy could probe certain classes of sources and/or 
test certain propagation models in presence of magnetic fields to be asso- 
ciated with such large scale celestial patterns. In addition, the evidence for 
large scale anisotropy around 1 EeV 1 claimed b y the AGASA collabora- 
tion ( Havashida et al. 1 998; H avashida et al. 19991 ) motivates even more this 



kind of studies. 

From the PeV energy range, the flux of cosmic rays is so low that we need 
ground based experiments with large collecting areas measuring the secondary 
products of the interaction of the cosmic ray in the upper atmosphere. As 
any ground based experiment has only at one's disposal a limited field of 
view in declination distribution, anisotropy analysis are generally done owing 
to the nearly uniform exposure in right ascension by using a 1-dimensional 
coordinate system instead of the natural 2-dimensional one over the sphere. 
This is the case for example of the Rayleigh formalism ( Linslev 19751 ) which 



necessarily corrupts the sensitivity to tiny anisotropies. In order to exploit the 
angular power spectrum analysis methods, it is assumed wi thin the cosmic 
rays community that a full exposure of the sky is required ( Sommers 2 000; 



Anchordoqui et al. 2003f ). The aim of this paper is to show that this conclusion 



arises only because of the choice of the spherical harmonic coefficients estimate, 
and to show that with another choice of estimate, standard anisotropy analysis 
methods can be used even with a partial and non-uniform coverage of the 
celestial sphere. 

By denoting Hi each cosmic ray arrival direction, the standard estimate of the 
spherical harmonic coefficients is computed through 

i A Y^ifk) 

at m = —r 2^ 



where N is the total number of events, uo{n) is the relative exposure func- 



1 1 EeV= 10 18 eV 



2 



tion of the considered experiment, and A a normalization constant taken as 
J2iLi l/<^(^j)- A s we U known, the use of 1/u allows for decoupling the modes 
when working with a variable exposure over the whole celestial sphere, but 
breaks down in case of partial exposure of the sky, because it is no longer 
possible to perform the full sky integrations tha t are required t o measure the 
multi-poles of the celestial cosmic ray intensity ( Sommers 2000[ ) . 



In this paper, we choose to introduce and adapt the quadratic estimator 
method that is widely use d in the Cosmic Microwave Background analysis 
(see e.g. ([Hivon et al. 20 02)) where effects of a partial exposure can be decon- 
voluted from the observations in order to recover the true underlying power 
spectrum. We show that the application of this method allows for the standard 
anisotropy analysis with an exposure possibly going to zero in some parts of 
the sky. This point is of major interest for most cosmic rays experiments, as 
the Southern site of the Auger Observatory for instance. 



Another app roach to power spectrum estimation is th rough maximum like 



lihood (see (|Bond et al. 19981 : iBorrill et al. 199 9) and ((Hamilton 20031 ) for a 



review) which has the advantage of solving exactly the problem. It however re- 
quires an explicit representation of the sky covariance and is computationally 
very time consuming and numerically hard to achieve on large datasets. The 
quadratic estimate proposed here avoids this difficulty using a Monte-Carlo 
simulation (this is why it is often refered to as a "frequentist approach"). 



This paper is organized as follows : in the next section, we describe our method 
and compute all the statistical properties of our choice of angular power spec- 
trum estimates. From section 3 and on we apply it to the forthcoming Auger 
observatory. In section 3, we present the relevant informations about this ex- 
periment that we need in the context of angular power spectrum estimation. 
In section 4, we discuss the constraints that the Auger Observatory can put 
on isotropic distribution of cosmic rays at ultra-high energy. At last, in section 
5, we extend the analysis to the case of a large scale anisotropic distribution. 



2 Angular power spectrum with a partial sky coverage 



2. 1 Generalities 



The number of cosmic rays observed per unit solid angle dN/dQ is a Poisson 
random variable in each direction n, whereas considered as a function of ft, this 
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is a Poisson random process. We model it with the two dimensional quantity 



dN 

dn 



N 



(n) = M{H) = ^ S(H, n,) 



where 5 is the Dirac delta function on the surface of the unit sphere, and Hi 
the position of the i th cosmic ray. The total number of cosmic rays observed 
is then / M{H)dH = N. This distribution follows a Poisson law V(u(H)) with 
an averaged intensity density in the direction H : 



where W is the relative coverage of the experiment varying from to 1, fi = 
J W(H)dH the fraction of the sky effectively covered by the experiment, and 
A some continuous stochastic field that measures the departure from isotropy. 
The stochastic field A is assumed to have a zero expectation value : 



where we have introduced the average over all the possible realizations of the 
random phases of the A field. The expansion of A on the spherical harmonics 
basis is given by : 



e>om=-e 

It is natural to introduce the A random field because of the probable stochas- 
tic nature of the cosmic ray sources distribution and propagation through 
magnetic fields. The possible anisotropics we want to characterize are a pri- 
ori of random nature. In particular, many theoretical models of UHECRs 
anisotropies are based on at least partly random configurations of sources and 
magnetic fields. We are driven to interpret a particular set of events as one 
specific realization of the random process; and we are led to characterize the 
underlying random process properties through the angular power spectrum of 
the data of the only set of events we have. 

Turning now to the two point correlation function of the A field, we assume 
this function to be only dependent on the angular distance between two points 
on the sphere : 



N 



W(H) (1 + A(n)) 



4tt/i 



(A(n)) r = 



(A(n)A*(n')) r = £(n • n) 
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This is a very strong and important hypothesis as it is the basis of any power 
spectrum estimation, in the classical sense we want to give to it. In particular, 
this assumes that the statistical properties of A are the same over the whole 
celestial sphere. In the following, we will refer to this hypothesis as spatial 
stationarity by analogy with a stationary time dependent problem, where the 
notion of stationarity is used when the two point correlation function depends 
only on the time difference 2 . Making a simple expansion of £(n • ft) onto the 
Legendre polynomials 3 : 

2/ 4- 1 

^(n ■ ft) = £ -^C e P e (n ■ ft) = £ C t Y tm {n)Y? m {n") 

£ £,m 

and of (A(n)A*(n')) r onto the spherical harmonics : 

(A(n)A*(n')) r = £ £ (a, imi <™ 2 ) r W^™ 2 (™') 

£i,mi £2,m,2 



we clearly see by identification that spatial stationarity leads to a diagonal 
covariance matrix of the a^ m coefficients : 



where Ce is the angular power spectrum of the fluctuations. 

The angular power spectrum is therefore a two point correlation function in i 
space. It gives information on the correlation between two angular directions 
separated by an angular scale ~ l/£ (in radians). For a Gaussian field A, 
the Ce power spectrum characterizes completely the fluctuations of the field 
as the even order moments are obtained from the second order moment and 
the odd order moments are zero (Wick's theorem). In the general case, this 
is no longer true and higher order moments (three points correlation function 
and so on ...) are necessary in order to fully characterize the field. One should 
note however, that, if the observed comic ray sky is the result of a sequence 
of many random processes and assuming no single process dominates (e.g. the 
sky is not the result of a single dominant source with no magnetic field), the 
fluctuations will be of Gaussian nature according to the central limit theorem. 

2 We use here the vocabulary of "spatial stationarity" rather than the one of "ho- 
mogeneity" by arbitrary choice. 

3 We use here the spherical harmonics addition theorem : 

£ Y em (n)Y; m (n>) = -^Pe(n ■ ft) 
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In any case, higher order moments can be computed on the data and departure 
from a Gaussian behavior can be measured. 

As outlined before, the data set we are dealing with is a Poisson sample of the 
random field A. Consequently, we have to introduce a second kind of average : 
the average over all possible sample configurations {-) p . Therefore, from now 
on, we use the notation (•) = ( (-) P ) to express this double average over the 
possible configurations of N events and over the possible realizations of the A 
random field. From the elementary Poisson statistic properties, it is easy to 
show that : 

(Af(n)) = (Mn)) P ) r = (u(n)} r = ^j-W(n) 



and : 



N n2 



(Af(n)Af(n')) = [^j W(n)W(n')(l + (A(n)A(n')) r ) 
AT 

+ W {n)8{n,n') 

47T/i 



2.2 Definition of the spherical harmonic coefficients estimate 



We want to build an estimate of the harmonic expansion coefficients of A. In 
the cases we are interested in, the field A is not measured uniformly over the 
whole celestial sphere. This is due to the non uniform exposure of cosmic ray 
experiments. For a single experiment, the knowledge is even limited to a given 
region in the sky and no information on A is available elsewhere. Moreover, in 
this given region, the exposure is not uniform and generally depends on dec- 
lination. When combining data from two observatories, the exposure becomes 
full sky but non uniform. This is show n for instance with Sugar and AGASA 
coverage in dAnchordoaui et al. 2003f ). or with Auger Southern and Northern 



sites in f Sommers 2 000). 



All these configurations can be described through the introduction of the win- 
dow field W(n) that measures the relative exposure in the direction n on the 
sky. This field can even vanish in some regions. Thus, A(n) = A(n) x Win) is 
the quantity we have access to experimentally and not simply A(n) as in the 
case of a uniform and full sky coverage. This has an immediate effect in the Ce 
determination as we cannot compute the expansion of the field we intended 
to. We only have access to what is called the pseudo-power spectrum Ce of 
the product of the two fields. A simple way to go back to the true Ce from the 
measurement of Cf was proposed for Cosmic Microwave Background analysis 
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by (IHivon et al. 2 002) and has been widely used in this community for variou s 



experiments (|Netterfield et al. 20021 iBenoit et al. 20021 iHinshaw et al. 20031 ) . 



We will soon show that the convolution kernel which mixes the modes of the 
angular spectrum we want to measure is the same as the one found in the 
framework of the CMB. 

We denote our estimates ag m and we define them as : 
a lm = J dnY im (n) 

4?r 47r /i 

Clearly, (a,£ m ) = 0, as well as for the expectation value of the true coefficients. 
2.3 The bias on the angular power spectrum estimate 



For reasons that will soon beco me clear, we introduce the following coupling 
kernel as in (|Hivon et al. 20021 ) : 



l\m\ 4^ 



K Mm > = V w hmi / dnY em ,(n)Y£ m (n)Y iimi (n) 



where we have expanded the window field on the spherical harmonics basis. 
Let us also introduce the power spectrum of the window field : 



m=—t 



fin 



Turning now to the correlation between two multi-pole estimates, it is easy to 
show that : 

i\m\ 



We then estimate the power spectrum Ci simply by taking the empiric average 
over m : 

■y m=£ 

Ct = 97TT ^ l^" 1 ' 

m=-l 
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This yields to : 



(: '' ' ~ 2f + 1 E f E \Kimi>imi\ 2 H — -Jf—Kemem 

m=—£ litni 



In ( Hivon et al. 2002^ . it has been shown that the first term is equivalent to 



a mode-mode coupling matrix Mu 1 : 

Ci) = E + 2i+ 1 ^ ^m«m 

£1 m=—£ 



where the M« x matrix elements are : 

2 



47r € 2 



1 ^2 





which makes use of the Wigner 3-j symbols. By expanding the second term 
onto the Wigner 3-j symbols, and after some manipulations, it is easy to show 
that : 

t i 



2 * + 1 m=-€ V47T 



leading to : 



C,) = ^M« 1 Q 1 + 



iV 



We therefore have a simple and analytical link between our estimate and the 
true Ci for a sky observed with a varying and/or incomplete exposure. Apart 
from a bias, our estimate is just the convolution of the true power spectrum 
by a kernel whose properties can be determined analytically from the shape 
of the window. 



At last, in ( Hivon et al. 20021 ). it is shown that the effect of M on a constant 



is a multiplication by the second moment of the window f2 — j~S W 2 (n)dn 
Ee ^W e . Therefore, we can go back to the angular power spectrum of the 
A field through : 

(cr p >=E^ 1 ( < ^) = ci+ w 
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We see that the experimental power spectrum is unmixed and asymptotically 
unbiased. The bias term can be easily computed analytically and is purely 
induced by the finite number of arrival directions that are available, that is, 
purely induced by the Poisson statistics of Af. 



2.4 The variance of the angular power spectrum estimate 



From the fourth moment of M and the Wick's theorem, there is no difficulty 
to compute the correlation between four multi-poles estimates. However, this 
calculation is rather long and tedious, so we don't reproduce it in details. As 
in cosmic ray physics, the null hypothesis we want to test is isotropy, we are 
interested in Cg = 0. In this case, the result for the covariance on Ct is found 
to be : 



N J 2£' + l 

Therefore, the variance on the experimental power spectrum simply reads : 



h,e 2 



2.5 Discussion 



One might ask the question of the pertinence of measuring an angular power 
spectrum on a partial region of the celestial sphere and ask what is the link 
between such a local power spectrum and the global one. In fact, all of this 
discussion is linked to spatial stationarity of the random field. If the field is 
spatially stationary, then, a partial part of the sky if a fair sample and can allow 
to recover all modes provided of course that the matrix M^t can be inverted, 
which is the case when the portion of the sky covered is larger than ~ 40% 
as shown on Fig. 1. If the sky is not spatially stationary (which can be seen 
by comparing power spectra in various partial regions of the sphere), then our 
procedure cannot be applied and the recovered power spectrum is of course not 
valid for the whole sphere. In that case anyway, the non stationarity of the sky 
makes the notion of a single power spectrum on the sphere totally irrelevant 
and our procedure allows one to construct power spectra for different regions 
provided the fact that they are spatially stationary enough. We therefore see 
that in any case, computing a local version of the angular power spectrum is 
interesting as it provides a global information when it is meaningful (spatially 
stationary sky) or a local one if relevant. 
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Fig. 1. Error bars of the reconstructed dipole as a function of covered portion of the 
sky. For an experiment covering more than 40% of the sky, the error bars become 
stable and the corresponding dipole estimate makes sense. For values smaller than 
40%, the M| 1 ^ 2 matrix is no longer regular and, as a result, the error bars explode. 
This conclusion arises with the conservative choice of # m ax=60 deg. See section 3 for 
full explanations about the parameter # max - Blue and red crosses show the error-bars 
for Auger South and Auger North respectively; the dotted the expected level for a 
full uniform sky coverage. 



There are a lot of theoretical motivations to detect large scale patterns around 
1 EeV as well as at higher energies. Magnetic fields or relative motion of the 
observer with respect to cosmic ray rest frame are natural mechanisms that 
lead to low-order moments. As these low-order poles characterize properties 
over the whole sphere, the condition of spatial stationarity should be naturally 
achieved in those cases, and local power spectra should show the same pat- 
terns whatever the observed portion of the sky. On the other hand, at higher 
poles, and if sources are located in large scale structure such as the Virgo 
cluster or the Super Galactic Plane, we expect that power spectra obtained 
from different regions of the sky covering or not those regions will show differ- 
ent characteristics. Again, spatial stationarity would be achieved within each 
region or along superstructures and the power spectrum obtained in those 
cases will give us informations on the sources distribution and magnetic fields 
within those structures. Of course, there is no way to completly describe the 
sky without looking at it. If a spectacular source or set of sources is present 
somewhere, the only way to know about it is to achieve (with several obser- 
vatories) full sky coverage. 
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The formalism of the angular power estimation is up to now the tool the most 
sensitive to search for tiny anisotropies over the sphere, provided the fact that 
the detector doesn't smear out the arrival directions more than the scale l/£. 
One can imagine that the use of this observable in the next future will be 
of great help in order to bring strong constraints on the UHECRs models of 
production and propagation. 



3 The Pierre Auger Observatory 



The Pierre Auger Observatory 4 is the first of a new generation of detectors 
specially dedicated to the highest-energy cosmic rays. Large area ground based 
detectors do not observe the incident cosmic rays directly but the Extensive 
Air Showers (EAS), a very large cascade of particles, that they generate in the 
atmosphere. All experiments aim to measure, as accurately as possible, the 
direction of the primary cosmic ray, its energy and its nature. There are two 
major techniques used. One is to build a ground array of sensors spread over 
a large area, to sample the EAS particle densities on the ground. The other 
consists in studying the longitudinal development of the EAS by detecting the 
fluorescence light emitted by the nitrogen molecules which are excited by the 
EAS secondaries. 



The Auger Observatory ( The Pierre Auger collaboration 1995f ) combines both 



techniques. The detectors are designed to be fully efficient for showers above 
10 EeV, with a duty-cycle of 100% for the ground array, and 10 to 15% for the 
fluorescence telescopes. The 1600 stations of the ground array are cylindrical 
Cerenkov tanks of 10 m 2 surface and 1.2 m height filled with filtered water; 
they are 1.5 km spaced on a triangular grid. 

The Auger observatory covers both hemisphere, one site in the South is 
presently under construction in Argentina and will be completed at the end of 
2005. The Northern site construction should start soon after. Once fully com- 
pleted in 2008, the Auger Observatory will be covering a surface of 2x3000 km 2 
and will provide unprecedented statistics. With a total aperture of more than 
14000 km 2 -sr, and with an integral cosmic ray intensity above 10 EeV of ap- 
proximately 0.5/(km 2 -sr-yr), the Auger Observatory should detect every year 
of the order of 7000 events above 10 EeV and 70 above 100 EeV (assuming 
a E~ z dependence of the spectrum). The angular resolution of the surface 
detector alone is believed to be of the order of 1 degree, and can be im- 
proved to less than 1 degree in the case of the so-called hybrid events, that 



Named after the French physicist Pierre Auger (1899-1993) who discovered the 
Extensive Air Showers. 
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is, e vents detected by both the surface d etector and the fluorescence detec- 
tor ( The Pierre Auger collaboration 1995f ). 



A full-time operation of the surface detector means that there is no exposure 
variation in sideral time and therefore constant exposure in right ascension. 
For such a detector located at latitude ao and fully efficient for cosmic rays 
arriving with zenith angles 9 less t han some maxim al value 9 msx , the exposure 
is only a function of declination 5 ( Sommers 20001 ) : 



W(n) oc cos (a ) cos (5) sin (a m ) + a m sin (a ) sin (<5) 



where a m is given by : 



o if e > i 

7T if f < -1 

cos" 1 (£) otherwise 



and : 

_ cos (6> max ) - sin (op) sin (5) 
cos (ao) cos (S) 



Fig. 2 shows the resulting declination dependence for the two sites of the 
Auger Observatory located at latitude ao = +39 deg. for the Northern one, 
and a = —35 deg. for the Southern one. The cut angle 9 max is chosen at 60 
deg. Also shown is the combined exposure, which will completely cover the 
celestial sphere but will not be uniform. 

In order to show that analysis are not sensitive to the choice of # max for the 
two latitudes we are considering, Fig. 3 plots, for a set of 7000 events, the 
variation of the first multi-pole C^ xp error bars as function of # max . Clearly, 
for both sites from 50 deg. and on, the error bars are stable with respect to 
the ^max parameter. 



4 Predicted constraints on isotropy with the Auger detector 



In this section, we choose to deal with events with energy beyond 10 EeV, 
where the required fully efficient cosmic rays detection is satisfied by the Auger 
Observatory for a large range of 9 ma _ x . The total number of events N tot (E) 
being detected with the Auger arrays on the whole sphere is approximately 
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Fig. 2. The declination dependence of the Auger Observatory relative exposures. 
The Southern and Northern sites are indicated separately by dots, whereas the 
combined exposure is in solid line. 
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Fig. 3. Amplitude of the C^ xp error bars as function of max for the two Auger sites 
latitudes, computed with a set of 7000 events. From # max = 50 degrees, results on 
the CT V error bars are not sensitive to the choice of # m ax- 
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N tot (E > lOEeV) = 7000 per year. We consider here the number of events 
N ev for an integration time T, weighted by the covered fraction of the sky : 

N ev = N tot (E> lOEeV) xTx f t 



To check our estimate of the bias and of the error bars, we simulated 100 times 
N ev events on a uniform (zero power spectrum) sky with the Southern Auger 
site coverage. We then reconstructed the power spectrum using the relations 
given in the previous sections. Fig. 4 shows the perfect agreement between the 
Monte-Carlo simulations and the analytical predictions for both the bias (the 
bias has been subtracted and the resulting values are indeed centered on zero) 
and the error bars. This plot is performed for an integration of 1 year of data 
taking. 
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Fig. 4. Comparison between the analytical estimation of the bias and of the er- 
ror-bars and the simulated ones of our angular power spectrum estimate. The cov- 
erage of the experiment is assumed to be the Auger Southern one, with a duration of 
1 year data taking. We consider only statistics beyond 10 EeV. Clearly, the analyti- 
cal computation perfectly reproduces the properties of the Monte-Carlo simulation. 

The enhancement of the number of events with the Northern site allows for 
more stringent constraints on an isotropic distribution as can be seen on Fig. 5. 
We have imposed 3 years of data taking for the Northern site, and 6 years for 
the Southern one (3 years for the only Southern site + 3 years for both the 
Southern and Northern sites). With such statistics, it becomes possible to test 
the isotropy hypothesis with an accurate precision. 
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Fig. 5. Constraints on an isotropic distribution of cosmic rays beyond 10 EeV for 3 
years of data taking of the Southern site only (dashed line) , and for 6 years of data 
taking of the Southern site + 3 years of the Northern one (continuous line). These 
error-bars are estimated from the analytical computation. 



5 Sensitivity of Auger to a dipole 



In order to check for the efficiency of our mode deconvolution, we simulated 
events following a uniform distribution plus a dipole of 10 % amplitude. In 
practice, from the Cg input spectrum, we generated the ae m coefficients with 
uniform random phases. We then transformed these coefficients into a sky map 
using an inverse harmonic transform in order to have a realization of the A 
random field. We then multiplied this map by the required coverage of the 
sky, and drew the number of events falling in each pixel using a Poisson law 
with average proportional to the this map. The exact position of each event 
within the pixel is drawn uniformly. All these steps rel y heavily on the soft ware 
provided along with the Healpix pixellisation scheme ([Gorski et al. 19 98). For 
each Monte-Carlo realization of a set of events for a given input power spec- 
trum, we then extract the pseudo power spectrum by expanding the map onto 
the spherical harmonics basis (using the anaf ast routine) and then apply the 
deconvolution to reconstruct the unbiased and unmixed power spectrum. We 
did the simulations for the Southern site only and for both sites with the same 
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integration of data taking than in the previous section. 



The result of the Monte-Carlo are shown on Fig. 6. The reconstructed power 
spectrum is in agreement with the input one, and excludes the isotropic distri- 
bution at the 2.5 a level using the Southern site only. Obviously, the improve- 
ment of the reconstructed power spectrum is clear with the 3 years data taking 
of the Northern site added to the analysis. As a consequence, the error-bars 
of the quadrupole are strongly reduced. The same procedure can be applied 
to any input power spectrum. 
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Fig. 6. Reconstructed power spectrum (expressed here in [Number of 
events/year/str] 2 )in case of a pure dipole input sky. The diamonds show the input 
power spectrum. In red is shown the reconstructed power spectrum with the only 
Southern site of the Auger Observatory, for a duration of 3 years of data taking; 
whereas in blue is shown the reconstructed power spectrum with both sites, for a 
duration of 6 years for the Southern site and 3 years for the Northern one (as in 
previous section). 



6 Conclusions 



We showed that in the general case of a varying and incomplete exposure on 
the sky, the true power spectrum of the cosmic ray sources distribution can be 
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recovered. This result is not new in itsel f as it was introduc ed a few years ago 
in the framework of CMB data analysis ( Hivon et al. 20021 ). Its application to 
cosmic ray data is however new and might open new possibilities as the general 
feeling up to now was that no power spectrum can be reconstructed with- 
out a complete sky coverage. The power spectrum that our procedure allows 
to recover is equivalent to the full sky one if the anisotropies in the arrival di- 
rections of the cosmic rays are well modeled by a spatially stationary random 
field on the sphere. If this is not the case, the recovered power spectrum is still 
valid, but only for the region that was used to determine it. Anyway in the 
non spatially stationary case, different power spectra are required in different 
regions of the sky and our approach is still relevant. Additionally, we have 
analytically solved the calculation of the bias and of the variance introduced 
by the finite sampling of the sky in the general case of a varying and eventu- 
ally incomplete exposure. Using the deconvolution proposed here, any cosmic 
ray dataset will be usable for anisotropy determination purpose, provided the 
fact that the arrival directions and coverage map are known within reasonable 
precision. 
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